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ABSTRACT 

We model the observed size and brightness of the VLBA radio core of the jet in Cygnus X-1 to derive 
an expression for the jet power as a function of basic jet parameters. We apply this expression to 
recent constraints on the jet power from observations of a large scale shocked shell around the source 
by Gallo et al. 2005, which leads us to a set of alternative conclusions: either (a) the jet contains 
large amounts of protons: (> 2000 protons per radio emitting electron), (b) it has a very low volume 
filling factor of / < 3 x 10~ 5 , (c) the steady, radio emitting VLBA jet is not the source of the kinetic 
energy powering the ISM shell, or (d) its asymptotic behavior differs fundamentally from a broad set 
of plausible analytic jet models. 

Subject headings: ISM: jets and outflows — X-rays: binaries — radio continuum: general — X-rays: 
individual (Cygnus X-1) 



1. INTRODUCTION 

The field of X-ray binary (XRB) study has seen a re- 
markable transformation over the past 5 years: Relativis- 
tic jets have moved from being considered exotic and 
rare abnormalities to being recognized as integral and 
maybe vital components in the transfer of energy and 
angular momentum in accreting stellar mass compact ob- 
jects. While it has been known for decades that the jets 
produced by supermassive black holes carry enormous 
amounts of energy and transform their environment, it 
has been difficult to make firm estimates of the kinetic 
power of XRB jets. 

Furthermore, the study of XRB jets is plagued by 
the same dilemma facing extragalactic jet research for 
decades: a lack of solid information about jet composi- 
tion (i.e., do jets mostly contain protons or positrons as 
positive charger carriers, and is the inertia dominated 
by non-thermal relativistic particles or cold, thermal 
plasma), jet speed (how relativistic are the jets in their 
very cores), and the jet power. These questions have 
been discussed for decades Ce.g. IRe ynoldse t all 119961 
IWardle et aflll99l IGhisellini fc Celottill200i]T "yet con- 
clusive observational evidence to answer them is sparse 
at best, so the debate continues. In this situation, tight 
constraints from individual objects are our best bet at 
making progress on any of these questions. 

The interest in XRB jets has arisen from the discovery 
of the nearly universal presence of compact, flat spectrum 
radio emission from accreting black hole X-ray binaries 
ijFenderl l2001|) in the so-called low/hard state, where 
most of th e X-ray emission is in the form of a hard pow- 
erlaw ( see iMcClintock fc Remillardl l2003. for a detailed 
review of black hole X-ray states). In two cases, this 
radio emission has been resolved into clearly collimated 
jets iphawan et al.H2000l: IsTirling et alJ200llfl . The spec- 
tral and morphological similarity to the cores of AGN 
jets further strengthens the interpretation of this radio 
emission as evidence for jets. 

However, these jets are not resolved transverse to their 
axes and in the absence of other information, it has been 



difficult to constrain the kinetic power of these jet - open- 
ing angles have been assumed, but the results are then 
always subject to an unknown, arbitrary parameter. Fur- 
thermore, the information encoded in the spatial distri- 
bution of the observed radio emission of the VLBI obser- 
vations of the two jets which have been resolved (Cygnus 
X-1 and GRS 1915+105) have not been harvested to a de- 
gree that allows tight constraints to be derived. The aim 
of this paper is to derive just such constraints. Cygnus 
X-1 is far better suited for such a study because its dis- 
tance is know to significantly higher accuracy than that 
of GRS 1915+105. 

The status of firm observational data on XRB jet power 
was also recently improved by the discover y of a shell 
of th ermal emission around Cygnus X-1 liGallo et alJ 
2005). This shell has been interpreted as the shocked 
ISM around a low surface brightness, jet-driven radio 
lobe. Using analysis borrowed from radio galaxy dy- 
namics, the authors limit the average kinetic power (W) 
from the source to fall between 3 x 10 36 ergs s _1 < (W) < 
3x 10 37 ergs s^ 1 . This discovery, together with the VLBA 
observations of resolved jet emission, makes Cygnus X-1 
an ideal candidate to tighten the constraints on funda- 
mental jet parameters. 

In order to put these limits into the context of the 
compact VLBA jet observed bv lStirling et alJ (|200lD . we 
apply a representative emission model to the VLBA data 
(a2J. This allows us to derive a parameterized expression 
for the kinetic power which depends only on a few un- 
knowns which we hope to constrain, namely the particle 
content of the jet, the filling factor of emitting material, 
and the equipartition fraction of the magne tic field. In 
jQwe compare the outcome to the results of l|Oallo et alJ 
2005), which is orders of magnitude larger than the es- 
timate we derive for our fiducial set of parameters, and 
discuss the quantitative constraints that can be derived 
on the set of interesting parameters from this compari- 
son. Section Q] summarizes our findings. 

2. THE EMISSION FROM THE JET IN CYGNUS X-1 
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The radio spectrum of the jet in Cygnus X-l is flat 
to slightly inverted (a = d\ogL v / d\ogv > 0), with a 
flux level varying betw ee n a few and a few te ns of mJy 
ijBrocksopp et allll999ft . iStirling et all 1)20011 epochs A 
and C) resolved about 50% of the emission to be ex- 
tended, measuring roughly \i ~ 10 — 15 mas. The canon- 
ical mo del for the flat radio emissio n from jet cores goes 
back to lBlandford fc Koenigll l)1979l BK hereafter). This 
model has been used successfully and extensively in the 
context of both AGN jet core s and X-ray binary je t s (e.g . 
iHiellming fc Johnston! Il988t iFalcke fc Biermannl Il996fl 
and we will employ it here as well. The underlying as- 
sumption of the model is that the jet is described well by 
a freely expanding flow with roughly uniform velocity, 
implying a constant half-opening angle <j). Such a de- 
scription should be a good first order approximation at 
least in a limited region around the location where most 
of the radio emission originates. In such a jet, the mag- 
netic field evolves quickly into an almost purely toroidal 
topology and follows B oc z , where z is the distance 
along the jet, measured from the jet origin at the com- 
pact object. 

Rather than assuming adiabatic evolution of the par- 
ticle distribution, the BK model imposes fractional 
equipartition between particles and magnetic fields such 
that the particle pressure is p pa rt = oc z~ 2 is a fixed 
fraction of the magnetic pressure p mag , where £b is the 
equipartition fraction between particles and the magnetic 
field 1 . One may allow for the presence of protons with a 
contribution pp r ot = £pPpart/(l +£ p ) to the particle pres- 
sure. The electrons (and possibly positrons, subsumed 
in the non-proton part of the particle pressure) are as- 
sumed to follow a powerlaw distribution with index s 
such that the number density follows 71(7) = C c j~ s . For 
simplicity and lack of better measurements in the case of 
Cygnus X-l we will use s = 2, which is the fiducial value 
often used in AGN jets. We will also assume that, due 
to Doppler boosting, the emission is entirely dominated 
by the approaching jet, which is supported by the lack 
of emission detected from the counterjet in Cygnus X- 
1. Finally, we assume that the synchrotron emitting and 
absorbing plasma has a volume filling fraction of / < 1. 

The synchrotron emissivity, measured in the observer's 
frame, is given by 



2.4 x 10~ 17 ergs2p 7 / 4 ^ /4 5 2 / 



cm 3 Hz s 
= Co piS 2 



1 + Ce 



3.4 GHz 



(1) 



and the self-absorption coefficient, again measured in the 
observer's frame, is 



2.3 x 10~ 12 2p 2 £ B 
cm 1 + £ p 
C lP 2 6 2 



8.4 GHz 



6 2 f 



with the obvious definitions 
l|Rvbicki fc Lightmanl 11979ft . p 

netic plus particle) pressure, measured in the jet frame. 



(2) 

of Co and C\ 
is the total (mag- 



1 £b is a free, unconstrained parameter. The rough proportion- 
ality between p par t and p mag near the radio emission region can be 
relaxed slightly. However, a strongly non-linear relation between 
the two would lead to steep synchrotron spectra and is thus ruled 
out observationally. 



5 = [r (1 — f3 cos (6))] is the Doppler factor 
velocity v 



The jet 



/3c = cy/l - 1/r 2 likely falls into the 
range of 0.5 < (3 < 0.7 and the viewing angle falls 
into the range 25° < < 50° iStirling et all l200ll 
iBrocksopp etaTl 12001 iGleissner et all 12004ft . We will 
use (3 ~ 0.6 and 8 ~ 35° as fiducial values, which gives 
5 = 2.5<5i.6- 

For radio frequencies, the innermost region of the jet is 
optically thick, not contributing significantly to the total 
luminosity (see Fig. 1), which allows us to extend the 
lower integration limit for the luminosity along the jet 
axis (denoted by the z-coordinate) to z — > 0. The lumi- 
nosity emitted in a frame comoving with the jet plasma 
is then 



L u = I dz sin (6) 
la 



ete^ (1 - exp (-tv)) (3) 



for the optical depth t u (z, x) = 



sin ( 



We 



then define the photospheric radius zq as the location 
where the optical depth through the spine of the jet is 
one, i.e., t v (zo,x = 0) = 1 and express all distances 
in units of zq and all other quantities relative to their 
value at zq. The pressure can then be written as p — 
Po{z/z )~ with 



po = Vsin(f?)/(2C 1( P<M (4) 
The jet luminosity from eq. © is then 

L„ = 5.1z ^ [sin (60] McbCrM (5) 

As was shown by BK, the emitted spectrum is flat, a = 0. 

The VLBA flux of the Cygnus X-l jet is F u ~ 12 mJy 
( Stirli ng et alll2001|) . For a distance of D = 2kpc-D 2 > 
the observed luminosity is 



i cyg nusx-i = 5.7 x lO^ergss^Hz- 1 !^ 2 



(6) 

Epochs A and C are very similar in their characteris- 
tics: About 50% of the flux is resolved along the jet and 
both have essentially equal total fluxes of 12 and 12.9 
mJy and equal angular resolution along the jet. Epoch 
B has a lower resolved- to-unresolved flux ra tio 

In epochs A and C of IStirling et~aTl l|2001ft . 50% of the 
jet flux is resolved along the jet axis with an angular 
resolution of /u 50 % ~ 3mas, we can go back to eq. (0 
and solve for the value of z^% for which 50% of the jet 
emission comes from regions with z > 250%. Numerical 
evaluation (see Fig. shows that this is the case for 
z 50% ~ 3zo. Given /J5o% ~ 3mas, we can thus determine 
the photospheric distance Zq: 



5.2 x 10 13 cm 



sin 35' 



sine 



M50% 



3mas/ 



)d 2 



(7) 



Epoch B (which has slightly lower resolution along the 
jet) shows a more compact jet, with only about 30% of 
the flux resolved, which would imply a somewhat lower 
value of zq (by about a factor of 2). This might be due to 
changes in viewing angle due to jet precession, or possible 
interaction with the powerful wind from the companion 
(see SJ3J). We will thus carry the dependence on (150% 
through for transparency. Strictly speaking, eq. Q is 
a lower limit on zq, since the VLBA observations might 
have resolved out and thus missed some of the larger 
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Fig. 1. — Cumulative luminosity profile L„(< z) from the core 
to z/zq (solid curve) and the optical depth to synchrotron-self- 
absorption across the jet at location 2/2:0 (dashed curve). Also 
shown is the location of the 50% radius, 2 50 % (dashed grey lines) 
at about z 50 % ~ 3zq. 



scale, low surface brightness emission. However, since 
the total flux in the VLBA image is close to the total 
flux seen by other instruments, it is unlikely that zq is 
off by more than a few tens of percent. 
It is noteworthy that, from eq. the flux density is 

F v = L u /(4nD 2 ) <x (m 50 %) ¥ D~* (an0)~^ (8) 

i.e., for a given angular size 1^50%, the flux from a com- 
pact jet is almost independent of the distance D to the 
source! Because the angular distance and the luminos- 
ity distance scale differently with redshift, it should be 
possible in principle to use jet cores as yard sticks. The 
main obstacle for such a use is the fact that relativistic 
beaming will introduce significant uncertainty and bias, 
such that a statistical analysis of this effect will almost 
certainly be flawed. However, given a distance estimate, 
a measurement of F v and ^50% will yield a constraint 

10 3 

on (sin 9) s 6~s . Together with independent constraints 
on 9 (e.g., from the jet-counterjet surface brightness ra- 
tio), this provides a constraint on the jet velocity, with 
very little error associated with distance uncertainties, 
which are dominant in veloc ity measurem ents based on 
superluminal proper motion (Fender 2003). 

By equating eqs. JSJ) and © and substituting from 
eqs. {2} and {JJ we can calculate the jet opening angle: 
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(9) 



The value we find for (ft is very small compared to jet 
opening angles observed in AGN jets (e.g., in the M87 
jet, the half-opening angle of the large scale jet is of order 
(ft ~ 1.5°), but entirely consistent with t he upper limit 
of (ft < 2° set by the VLBA observations iStirling et all 
(2001). If the jet is in free expansion, this small opening 
angle would imply a very low internal sound speed, i.e., 
a pool of cold particles in addition to the synchrotron 
emitting electrons (see 331 . Note that the value of (ft is 
rather sensitive to ^50% > which introduces significant un- 
certainty given the variance of this value between epochs 



A through C of St irling et alJ l)2001[) . up to a factor of 5 
(the largest value of (ft ~ 0.3° is found for epoch B). 
Inserting this into eq. (QJ we get 



Po' 



1 erg 1 + £ B 
cm 3 ofit fi 



sing 1+£ P M50% 12mJy 
sin 35° 2/Z?2 3 mas Fg.4GHz 



(10) 



Given the opening angle and pressure of the jet, we 
can now estimate the kinetic jet power for the leptonic 
and electromagnetic jet content (i.e., not counting the 
rest mass kinetic energy in possible thermal particles the 
jet might be carrying, but including both jet and coun- 
terjet): 



W min ~2 xlpY 2 (3cK{<f>zf 



•1.9 x lO^ergss^ 1 



i + £b 



r i.25^0.i 
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sin(0) / 12mJy 3mas\ 2 £b (1 + £p) 
sin (35°) ^F 8 .4ghz /i50%/ ~VDM7 

This estimate is significan tly smaller than the W ~ 
10 35 ergss _1 estimated by ISoencer et al . ( 2001) 2 and 
typic al values assumed in the literature for compact XRB 
jets l)Fender et al.l I2003T) . This is mainly due to the 
smaller opening angle required by the observed value of 
/i 50 % . Unlike the value for (ft from eq. the W estimate 
in eq. I|ll|l is less sensitive to the value of /iso%- The 

/Zgpy o dependence still introduces an uncertainty of up to 
a factor of 3, given the variance in uf,ci% and Fraghz 
seen b etween the three different epochs of Stirling ct al. 
l)2001|) . with epoch B leading to a power estimate of 
W min ~ 6 x 10 33 ergss _1 . 

3. DISCUSSION 

At face value, this result implies that the kinetic 
jet power carried by fields and radio emitting elec- 
trons is smaller than the radiative power in X-rays by 
about 4 orders of magnitude. If other X-ray binary 
jets in the low-hard state were similar in nature to 
the jet in Cygnus X-1 , this would imply that jets do 
not carry away most of the accretion power in low- 
ha rd state and quiesce nt sources, as recently suggested 
bv iFender etaP {2003). This would require an ineffi- 
cient accretion scenario like an ADAF {Naravan fc Yil 
11991. CDAF lIQuataert fc Gruzinovll2000D . or ADIOS 
( Bla ndford fc Beeelmanl I1999J) . It would also imply 
that recent, independent estimate s of the jet power 
from steady low/har d state jets ijFender et alJ 120051 
iHeinz fc Grimrr]l2005j) are off by four orders of magni- 
tude. 

Such a low power outflow in Cygnus X-1 is in direct 
contradiction to the energy requirement for the recently 
found therm al shell around th e putative radio lobe of 
Cygnus X-1 {Gallo ct aLl l2005j) . which requires an aver- 
age jet power of 

3 x 10 36 ergs s" 1 < W sheU < 3 x 10 37 ergs s" 1 (12) 

Furthermore, several low/hard state XRBs show IR 
excesses that have been convincingly interpreted as 



2 semi-analytic models by Markoff ct al. 1 2001 
icantly larger values for the power. 



also give signif- 
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jet emi ssion (XTE J 1 118+ 4 80. IFender et all 120011 GX 
339-4, iHoman et all pOOl iNowak et alJ 120051 GR.S 
1915+105. lQglev et al.ll2000]) with total radiative energy 
output from the jet in excess of 5% of the total radia- 
tive output, amountin g to L\ Pt > 2 x 10 3 5 ergss~ 1 in the 
case XTE .11118+480 l)Fender et al.l200lD . The jet power 
must exceed the radiative power by a significant margin 
since otherwise the jet would radiate away more energy 
than it is carrying and thus violate energy conservation 3 . 
While the bright companion of Cygnus X-l makes direct 
observ ations of an IR/opti c al con tribution by the jet dif- 
ficult, IFender fc Kuulkersl {2001) argues that a flat jet 
spectrum extending to the V-band would still require a 
kinetic power in excess of ~ 10 35 ergss _1 . 

It is thus rather unlikely that the power of the Cygnus 
X-l jet is truly as low as implied by eq. 1)11)1 for the fidu- 
cial choice of parameters. We will discuss four alterna- 
tive (not necessarily mutually exclusive) conclusions that 
can be reached from the mismatch between eqs. (jl 1|> and 

3.1. Jet composition 

It is clear from eq. (tTT|) that W m i n depends on £b and 
£ p , thus, jet composition and field strength both factor 
into the total power. 

It is straight forward to show that the a significant 
deviation from equipartition between particles and mag- 
netic field will not change the estimate of eq. (|11|) any- 
where close to what is needed to bring it in line with 
the limits from lGallo et al.l l|2005D . The asymptotic de- 
pendence on £b is as follows: for small £b, Wmi„ actu- 
ally decreases rather than increasing, thus exacerbating 
the power discrepancy. For large £b on the other hand, 
Wmin oc £g^ 9 , requ iring £b ~ 10 36 to bring W m i n in line 
with the results bv lGallo et alJ {2005). Such large values 
are clearly unreasonable and unphysical. In other words: 
we cannot constrain the magnetic field strength, but sig- 
nificant deviations from equipartition will also not bring 
eq. (111(1 in line with eq. I)12|l . 

On the other hand, if the jet contains a large amount 
of cold, thermal plasma, the kinetic power can be sig- 
nificantly enhanced. Before going into more detail, it 
is important to note that simple charge balance of cold 
protons with the observed synchrotron emitting electrons 
will not enhance the power by more than two orders of 
magnitude because the above estimate is based on mini- 
mum energy arguments. For £b ~ 1, the Lorentz factor 
of the electrons contributing most of the 8.4 GHz flux 
is about 7c.8.4GHz ~ 70 (po/lergcm -3 ), thus adding one 
cold proton per emitting electron would only raise the ki- 
netic power by about an order of magnitude, insufficient 
for bringing it in line with eq. 1121) . 

The kinetic jet power including the inertial term is 

W = W min + 2pc 3 T/3 (r - 1) 7T (<pz) 2 (13) 

where w — 7/(7 — l)p + pc 2 is the enthalpy. The amount 
of cold gas necessary to bring the power estimate in line 
with the large scale constraints from eq. 1)121 parame- 
terized as Wtot = 3 x 10 36 ergs/sW 36 .5, bv lOallo et all 

3 Even the X-ray emissio n in the low/hard sate systems might 
parti ally originate in the jet IMarkoff et alJBOOlHMarkoff fc Nowald 
2004) , in which case the power requirement would increased by an- 
other 2 orders of magnitude, though synchro tron emissio n is prob- 
ably not the dominant source of the X-rays IHeindl2004T) . 
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3/16 



r/3(r - 1) 



(14) 



or about 2000 protons per radio emitting electron (we 
have abbreviated the dependence on the underlying pa- 
rameters £bj £pj and / in the opening angle - see eq. ^ 
for the detailed parametric dependence). 

This is a firm lower limit. We should note again that 
the question whether jets contain large amounts of pro- 
tons has been discussed in the literature many times over 
but remains essentially unanswered to this day, espe- 
cially in the newer and less well studied class of XRB 
jets. Any new constraint on an individual object, like 
the one derived here, will thus be valuable addition to 
the debate. Note also that energetic requirements for 
the presence of cold, thermal protons (derived here) are 
different from the requirement of a thermal parent pop- 
ulation of leptons at the base of the jet, as required in 
models that attempt to match not only the radio but 
also IR through X-ray observations of X-ray binaries, 
such as those bv Iparkoff et al.ll20?)ll IMarkoff fc Nowakl 
120041 lYiian fc Cuill2005lL These particles are required on 
spectral grounds, not for energetic reasons and no quan- 
titative requirement for the presence of cold protons can 
be made from them. 

Arguments for the presence of cold protons have been 
made on dynamical grounds in order to explain the open- 
ing angles of je ts in the case of freely expanding, ballis- 
tic je t models (|Falcke fc Biermannlll996i IMarkoff et alJ 
1200 However, the opening angle and the jet power in 
those models are free parameters and thus constraints on 
the presence of protons in the jets are qualitative only. 
We will employ a similar argument in the following para- 
graphs, based on the new constraint we derived for the 
jet opening angle in eq. (jHJ. 

Above and beyond the lower limit from eq. (fHjl . the to- 
tal amount of cold plasma that might be traveling down 
the jet cannot be constrained in the absence of direct 
radiative signatures, and thus the power enhancement 
that could be provided by cold particles is formally un- 
constrained. However, the small opening angle implied 
by eq. © could be interpreted as evidence for the pres- 
ence of large amounts of cold matter if we assume that 
the jet opening angle is approximately equal to the Mach 
cone of the jet 4 , i.e., 



~ c s /{(3cT) ~ v/5p/3p/(/3cr) 



(15) 



Thus, if cold, thermal particles, such as protons (with the 
appropriate number density of cold electrons for charge 
balance), are present in the jet and are the reason for the 
small opening angle, we can use the estimate of <j> and 
eq. 115|) to obtain an estimate of the rest mass density p 
inside the jet. This is, in effect, an upper limit because 

4 Note that in a BK type jet the "temperature" (p/p) of the rel- 
ativistic plasma is constant along the jet, thus defining a tempera- 
ture based on the jet opening angle is well posed and meaningful. 
However, even if the thermal plasma were not of constant tempera- 
ture along the jet because of radiative cooling, this argument would 
still hold, as the opening angle would define the temperature at the 
location where the jet becomes ballistic (i.e., inertia dominated). 
Further out, the thermal temperature might be lower than at that 
point, but the opening angle would remain at its ballistic value. 
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the jet might be partially collimated by magnetic fields 
and thus not in free, ballistic expansion. 
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3.2 x l(T 15 g 
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1 ergs cm 



0.06° 0.75 

4> r/3 



(16) 



The jet power is then 



r/?(r- 



W<3 x lO^ergss" 1 — — ^ / P ,' (17) 

~ to 4.5 x 10- 16 gcm- 3 3/16 v ; 

In this case, the high particle densities of n ~ 2 x 
10 9 cm~ 3 might imply a detectable thermal X-ray flux 
from these particles (as is the case in the jet of SS433). 
The temperature of this thermal gas would be T = 
2.3 x lO 6 K(0/O.O6°) 2 (/3r/O.75). The total thermal 0.5- 
10 keV luminosity from the jet would be Lbrems ~ 
8 x 10 37 ergss _1 , comparable to the bolometric luminos- 
ity believed to be coming from the accretion flow. The 
associated soft X-ray/UV/optical line emission, which 
would be very similar to the spectrum of SS433, has not 
been observed. 

Such a dense plasma would cool rapidly 5 if it were 
present from the very base of the jet (the cooling time 
there would be about 10~ 9 seconds, much shorter than 
the dynamical time). In order to avoid catastrophic ra- 
diative cooling, particles would have to injected at a large 
distance from the black hole, z coo \ > 10 13 cm (where the 
radiative cooling time is longer than the dynamical time) . 
Coincidentally, this is comparable to the orbital separa- 
tion of ~ 10 12 cm of the system. 

It would be natural, then, to interpret the presence 
of such a putative thermal component as mass loading 
from the powerful wind blown by the companion. Since 
this thermal plasma is supposed to set the opening angle 
of the jet, collimation and acceleration would still have 
to be happening at the location where the mass load- 
ing takes place, which is very far from the black hole. 
This would provide interesting asymptotic constraints for 
models of jet dynamics. 

3.2. Low filling factor 

Another way to increase the power estimate from 
eq. I|ll|) would be a filling factor / much smaller than 
one. While the opening angle <fi depends only weakly 
on /, the jet power is somewhat more sensitive. A very 
low filling factor would translate to a physical picture 
where small regions dissipate a large amount of energy 
into relativistic particles, such as would occur in mag- 
netic reconnection events in a turbulent flow. 

In this context, it is interesting to note that high res- 
olution radio images of AGN jets (which have are much 
better suited for high resolution imaging due to their 
larger ratio of object size to distance) show very non- 
uniform and sometimes unresolved, filamentary emission 

5 It should be noted that a fixed opening angle in a cooling jet 
is not a contradiction: the opening angle of a ballistic jet is set 
at the point where the jet goes out of transverse causal contact, 
i.e., where the sideways expansion becomes faster than the effective 
internal sound speed of the jet plasma (e.g., the fast magnetosonic 
speed in a magnetized jet). Beyond that point, the jet will expand 
at an essentially fixed rate, possibly cooling adiabatically and/or 
radiatively. In other words, the temperature does not need to be 
proportional to <j> 2 everywhere along the jet. Thus, contrary to a 
first glance educated guess, no mechanism is necessary to reheat 
the plasma to maintain a fixed opening angle. 



at least in the regio ns where the jet is optically thin (e.g. 
IBiretta et aU ll995). It might well be that the emission 
from the smaller, optically thin region is similarly fila- 
mentary. 

In order to bring the esti mate of eq. Ill 111 b ack in line 
with the macroscopic limits (Ga llo et al.l2005l) . the filling 
factor would have to be / < 3 x 10 -5 , which is very 
small inde ed and difficult to achieve with internal shock 
scenarios llYuan fc Cuil^QQ Sj) or mu lti-zone continuous 
emission models ijMarkoff et aUfeOOlj) . 

3.3. A non-radiative source for the kinetic power 

It is also possible that the flat spectrum, steady, 
low/hard state jet that is observed directly in the VLBA 
radio images is not the source of the power that drives 
the ISM shell. 

For example, it is possible that the radio emission origi- 
nates in a spine along the jet axis, which is embedded in a 
much wider sheath of low emissivity plasma with opening 
angle s hccth- In this case, the true kinetic power could 
be much larger (by a factor of (^sheeth/^) 2 - To bring the 
power estimate in line with the value from iGallo et alJ 
( 2005), we would require </> s hccth > 2°. Note that this is 
different from a low filling factor scenario since the opti- 
cal depth through the spine is larger than that through 
a low filling factor jet of equal emission measure. Why 
most of the flow would be non-radiative in this case is 
similarly difficult to understand as the case of a low fill- 
ing factor, however. It would amount to postulating an 
unrelated, much larger angle flow that would supply all 
the power but no radiation and, in essence, would mean 
that the actually observed jet has no physical relation 
to the large scale bubble observed around Cygnus X-1. 
Whether a radiative signature of such a larger angle out- 
flow exists and how it would tie in with ADIOS and disk 
wind models is not clear. 

Since such a putative outflow has so far not been ob- 
served, there are no constraints on the nature of the 
flow other than the kinetic power one would have to 
postulate based on the large scale requirements. It 
would thus be impossible to derive mass flow rates and 
outflow velocities. However, given the observed power 
and using the Eddington accretion rate of M — 1.5 x 
10 19 g/s(M/10M Q ) for a 10 solar mass black hole as an 
upper limit on the mass flow rate, we can conservatively 

put a lower limit of v s hccth > 6000 km s~ l W^g 5 on the 
outflow velocity. 

It is also possible that most of the kinetic power in the 
Cygnus X-1 jets is supplied by the type of optically thin 
radio outbursts observed in many Galactic XRB sources 
like for example, GRS 1 915+105 l|Mirabel fc Rodriguez! 
ll994tlFe"nder et al.ll999fl rather than the steady compact 
jets discussed here. This would have interesting conse- 
quences for the estimates of duty cycles and radiative 
efficiencies for such events. Since they are unrelated to 
the VLBA observations discussed here, we will not dis- 
cuss this possibility in more detail. 

It should be noted that the flat spectrum IR emission 
from low/hard state XRBs that has been interpreted suc- 
cessfully as jet emission still places power constraints on 
the jet that, while smaller than eq. Ijl2|l . still far exceed 
that of eq. 111|) . IR observations of Cyg X-1 are difficult 
because of the high Galactic background and the bright 
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companion, thus we cannot apply these limits directly. 
However, if the Cyg X-l jet exhibits IR fluxes similar to 
those from other low/hard state jet sources, an unrelated 
source of kinetic power that might be responsible for the 
large scale ISM shell could not explain the tight IR/radio 
correlation, and thus the constraints derived in HH.ll and 
13. 21 would hold. Ongoing Spitzer IR monitoring will help 
to answer this question. 

3.4. Grossly different kinematics 

Finally, the treatment used in this paper might not be 
applicable because some of the basic assumptions we em- 
ployed might be false. This would imply that the most 
commonly used description for jet cores that captures 
the critical observable features, the "Blandford/Koenigl" 
model, would grossly misrepresent the kinematics and 
emissivity characteristics of jet cores (by many orders of 
magnitude). This in itself would be an important con- 
straint on jet modeling. We shall briefly consider whether 
a : 

The most obvious deviations from a BK model would 
be (a) a non-conical jet, (b) a deviation from the as- 
sumed p oc z~ 2 scaling, or (c) an electron spectral index 
s different from s = 2. Barring seriously pathological 
jet geometries and kinematics, the jet should be well de- 
scribed (at least near the radio photosphere) by power- 
laws for the jet radius R = Rq{z/zq)^ (with < C < 1 
for collimation) and the pressure p = po(z/zo) x . Fol- 
lowing eq. I|17fl . the kinetic power of the jet is W = 
[4pT 2 f3c + p (r 2 - T) /3c 3 ] ttR 2 . If the jet is not cold mat- 
ter dominated (pc 2 <C p) , the energy flux is only constant 
along the jet if x = ~ 2£, which we impose as a condition 
in the following. If it is cold matter dominated, once 
again i !3. II applies. 

The spectral index for such a jet would be a = 
(18£ — 8 — 2s — 3sQ/12(. Since the observed spectrum 
in low/hard state XRB jets is flat to slightly inverted (in 
the case of GX 339-4, a ~ 0.1, while at low frequencies, 
the jet in XTE J1118+480 is rather steep, a < 0.5), we 
can impose a > and find s < (18£ - 8)/(2 + 3£) < 2. 
Fermi acceleration, which is often invoked as the origin of 
the electron powerlaw spectra, typically produces spectra 
with s > 2, which is consistent with the above condition 
only for the case s — 2 and £ = 1, i.e., a BK jet. Slightly 
smaller values of s might be possible to achieve wit h rela- 
tivistic Fermi acceleration ijKirk fc Heavensj [l989') or re- 
connection, but the parameter space is severely limited 
to a region around the classical BK model (s = 2 and 
( = 1). Estimates presented in this paper will not be 
affected too strongly by such a small deviation in pa- 
rameters. 

One might also argue that the assumption of uniform 
velocity is incorrect, i.e., T might be increasing signif- 
icantly inside the radio emission region through active 
hydrodynamic or MHD driving, in which case a sim- 
ple kinematic model would be inappropriate. It is be- 
yond the scope of this paper to present dynamical jet 
models in any detail. However, it is clear that the ra- 
dio emission region is very far away from the black hole, 
roughly 10 7 gravitational radii according to eq. Q. Jet 
acceleration at such distances would be very surprising 
indeed, presenting serious asymptotic challenges to dy- 
namical jet models. Significant acceleration over such 



large scales would also imply enormous terminal Lorentz 
factors which would be at odds wit h the observational ve - 
locity constraints on Cygnus X-l IGleissner et al.l l2004). 

Thus, it is not trivial to construct a jet model that 
produces a flat spectrum, obeys energy conservation, has 
reasonable terminal velocities, and still deviates signifi- 
cantly from a BK-jet. 

However, it is worth keeping in mind that, while the 
Cygnus X-l jet is generally well described by a BK jet, 
it does exhibit some features that go beyond this sim- 
ple, robust pi cture: The jet does sh ow signs of curvature 
in epoch A of iStirling et al.1 l)200H) . which might be due 
to precession induced by the binary orbit, possibly ex- 
plaining the different appearance of epoch B compared to 
epochs A and C. Furthermore, other Galactic jet sources 
show an even broader variety of behaviors, such as spec- 
tral pivoting ijCorbel et all200"o|) and short term variabil- 
ity. Such behavior clearly requires more complex models, 
including the possibly confining and contaminating pres- 
ence of stellar winds (which, interestingly, should affect 
the Cygnus X-l jet most strongly, given the powerful 
wind from it 09.7 lab companion). Finally, the question 
of what accelerates the relativistic electrons inside the 
jet (as postulated by all jet emission models), which di- 
rectly relates to the assumed quasi-equipartition between 
particles and magnetic field in the BK model, remains 
unanswered. 

It is beyond the scope of this paper to present mod- 
els addressing these issues, but it is possible that more 
complex dynamical models, incorporating sophisticated 
particle transport calculations, will lead to a significantly 
different power estimate that is more in line with the ob- 
served values from eq. (fl2^l . As it stands, the require- 
ment to satisfy the observed flat spectrum and the VLBI 
surface brightness distribution (in detail), while at the 
same time satisfying the large scale power requirements, 
presents a benchmark for new jet models to measure up 
to. 

4. SUMMARY 

We presented a detailed analytic model of the VLBA 
emission from the Cygnus X-l radio jet to derive a para- 
metric expression of the jet power Wj c t- By applying this 
expression to new, independen t observational constraints 
on Wjet bv lGalloet al. ( 200^), we derived stringent con- 
straints on fundamental parameters of the jet, which lead 
us to the following set of alternative conclusions: The jet 
must either (a) contain large amounts of cold protons 
and/or (b) have am extremely low filling factor, and/or 
(c) the power that drives the ISM shell is not carried by 
the compact VLBA radio jet but some other, unrelated 
type of outflow (e.g., optically thin radio flares or a broad 
accretion disk wind) and/or (d) a strong asymptotic de- 
viation from the analytic Blandford-Koenigl model we 
employed (however, still subject to the very restrictive 
spectral and imaging constraints provided by the data). 
Each of these alternatives presents valuable input for jet 
modeling efforts. 
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Rob Fender for helpful discussions. Support for this work 
was provided by the National Aeronautics and Space Ad- 
ministration through Chandra Postdoctoral Fellowship 
Award Number PF3-40026 issued by the Chandra X-ray 



The Cygnus x-1 jet 



7 



Observatory Center, which is operated by the Smithso- National Aeronautics Space Administration under con- 
nian Astrophysical Observatory for and on behalf of the tract NAS8-39073. 



REFERENCES 



Biretta, J. A., Zhou, F., & Owen, F. N. 1995, ApJ, 447, 582 
Blandford, R. D. & Begelman, M. C. 1999, MNRAS, 303, LI 
Blandford, R. D. & Koenigl, A. 1979, ApJ, 232, 34 
Brocksopp, C, Fender, R. P., Larionov, V., Lyuty, V. M., Tarasov, 

A. E., Pooley, G. G., Paciesas, W. S., & Roche, P. 1999, MNRAS, 

309, 1063 

Brocksopp, C, Fender, R. P., & Pooley, G. G. 2002, MNRAS, 336, 
699 

Corbel, S., Fender, R. P., Tzioumis, A. K., Nowak, M., Mclntyre, 

V., Durouchoux, P., & Sood, R. 2000, A&A, 359, 251 
Dhawan, V., Mirabel, I. F., & Rodriguez, L. F. 2000, ApJ, 543, 373 
Falcke, H. & Biermann, P. L. 1996, A&A, 308, 321 
Fender, R., Maccarone, T., & Kesteren, Z. 2005, MNRAS, in press, 
press 

Fender, R. P. 2001, MNRAS, 322, 31 
— . 2003, MNRAS, 340, 1353 

Fender, R. P., Gallo, E., & Jonker, P. G. 2003, MNRAS, 343, L99 
Fender, R. P., Garrington, S. T., McKay, D. J., Muxlow, T. W. B., 

Pooley, G. G., Spencer, R. E., Stirling, A. M., & Waltman, E. B. 

1999, MNRAS, 304, 865 
Fender, R. P., Hjellming, R. M., Tilanus, R. P. J., Pooley, G. G., 

Deane, J. R., Ogley, R. N., & Spencer, R. E. 2001, MNRAS, 322, 

L23 

Fender, R. P. & Kuulkers, E. 2001, MNRAS, 324, 923 

Gallo, E., Fender, R., Kaiser, C, Russell, D., Morganti, R., 

Oosterloo, T., & Heinz, S. 2005, Nature, 436, 819 
Ghisellini, G. & Celotti, A. 2001, MNRAS, 327, 739 
Gleissner, T., Wilms, J., Pooley, G. G., Nowak, M. A., Pottschmidt, 

K., Markoff, S., Heinz, S., Klein-Wolt, M., Fender, R. P., & 

Staubert, R. 2004, A&A, 425, 1061 



Heinz, S. 2004, MNRAS, 355, 835 

Heinz, S. & Grimm, H. 2005, ApJ, accepted 

Hjellming, R. M. & Johnston, K. J. 1988, ApJ, 328, 600 

Homan, J., Buxton, M., Markoff, S., Bailyn, C. D., Nespoli, E., & 

Bclloni, T. 2005, ApJ, 624, 295 
Kirk, J. G. & Heavens, A. F. 1989, MNRAS, 239, 995 
Markoff, S., Falcke, H., & Fender, R. 2001, A&A, 372, L25 
Markoff, S. & Nowak, M. A. 2004, ApJ, 609, 972 
McClintock, J. E. & Remillard, R. A. 2003, ArXiv Astrophysics 

e-prints 

Mirabel, I. F. & Rodriguez, L. F. 1994, Nature, 371, 46 
Narayan, R. & Yi, I. 1994, ApJ, 428, L13 

Nowak, M. A., Wilms, J., Heinz, S., Pooley, G., Pottschmidt, K., 

& Corbel, S. 2005, ApJ, 626, 1006 
Ogley, R. N., Bell Burnell, S. J., Fender, R. P., Pooley, G. G., & 

Waltman, E. B. 2000, MNRAS, 317, 158 
Quatacrt, E. & Gruzinov, A. 2000, ApJ, 539, 809 
Reynolds, C. S., Fabian, A. C, Celotti, A., & Rees, M. J. 1996, 

MNRAS, 283, 873 
Rybicki, G. B. & Lightman, A. P. 1979 (New York: Wiley) 
Spencer, R., de la Force, C, Stirling, A., Garrett, M., Fender, R., 

& Ogley, R. 2001, Ap&SS, 276, 255 
Stirling, A. M., Spencer, R. E., de la Force, C. J., Garrett, M. A., 

Fender, R. P., & Ogley, R. N. 2001, MNRAS, 327, 1273 
Wardle, J. F. C, Homan, D. C, Ojha, R., & Roberts, D. H. 1998, 

Nature, 395, 457 
Yuan, F. & Cui, W. 2005, ApJ, submitted 



